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By applying a scattering-wave theory, the electromagnetic response of an ar- 
bitrary array of multiple slits perforated on a metallic film and filled with 
different slit dielectric materials can be studied in an analytical way. Here, the 
wavelength-dependent splitting of a light beam into two by asymmetrically 
filled slits in a metal film using intra- and inter-slit dual-wave interferences is 
fully explored. We consider a triple-sht structure perforated on a gold film, 
where the middle slit is used for the surface-plasmon excitation by a narrow 
Gaussian beam while the two side slits arc used for the detection of a trans- 
mitted surface-plasmon wave propagated from the middle opaque slit either 
at a particular wavelength or at double that wavelength, respectively. For this 
proposed simple structure, we show that only one of the two side observa- 
tion slits can be in a passing state for a particular wavelength, but the other 
blocked sht will change to a passing state at double that wavelength with a 
specific design for the slit depth, slit dielectric, and inter-slit distance in the 
deep sub-wavelength regime. In this sense, surface-plasmon mediated light 
transmission becomes wavelength sensitive in our model, and a single light 
beam can be separated into two according to its wavelength in the transverse 
direction parallel to the array. This provides us with a unique way for direct 
optical reading in the near-field region using a non-spectroscopic approach. 
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1. Introduction 

Surface-plasmon-polariton modes ^ and localized surface plasmons'^', which are both local- 
ized at an interface between a bulk conductor and a bulk dielectric, have become very hot 
research subjects in recent years ^. The extraordinarily high transmission of a p-polarized 
light beam propagating through a two-dimensional honeycomb lattice of holes on a metal 
film with sub-wavelength diameters'^ depends strongly on the lattice constant and the 
metal-film thickness (also in the deep sub-wavelength regime)'^. 

On the other hand, studies on the surface-plasmon mediated light transmission by a 
sub-wavelength structure (including random surface roughness) on a designed metal surface 
have also attracted a lot of attention'^nS'. For an optically-opaque metal film on a dielectric 
substrate, if a single slit is perforated on this film, the excited surface plasmons on the front 
side of the film can be coupled to the backside by intra-slit interferences'^'^. In addition, 
we demonstrate in this paper that not only the intra-slit interference but also the inter- 
slit interference for a slit array can affect the transmission of the excited surface plasmon 
propagating through the slits, which was not explored in previous worksl^l^. The inter- 
slit interference considered in this paper is related only to the surface wave but not to 
the surface-plasmon wave, since this interference effect survives even for a perfect electric 
conductor. It is also important to mention that the surface-plasmon-polariton mode for a 
planar surface becomes folded with a finite lattice constant in a periodic array and is split 
into many branches with a minigap opened either at the center or at the edge of the first 
Brillouin zone'^. Therefore, the surface-plasmon mediated splitting of a light beam studied 
in this paper has to consider avoiding these minigaps. A related work on plasmonic photon 
sorters for spectral and polarimetric imaging was reported early'^, which employed coupled 
bull-eye structures with a linear modulation in groove depth within each structure. 

In our previous study"^', we have shown the longitudinally color-dependent light focusing 
by a finite linear array of grooves with various widths in a parabolic pattern, where various 
focal lengths of a slit-array aperture were obtained for an incident plane wave with different 
colors. Here, rather than using a Green's function formalism^^HSl^ -^g present a scattering- 
wave theory which utilizes a slit-eigenmode expansion to treat an arbitrary array of slits 
having arbitrary spacings, widths and dielectrics. In addition, the derived scattering-wave 
theory in this paper provides one with a full description to the surface scattering of light by 
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removing a so-called "diagonal" approximation adopted in a previous related theory'^'^^J^. 
Specifically, as an example, in this paper we consider a triple-slit structure in which the 
middle opaque slit is used for the front-side local surface plasmon excitation by an incident 
Gaussian beam either at a particular wavelength or at double that wavelength, and the two 
side slits are used for the wavelength-dependent detection of surface-plasmon mediated light- 
beam splitting in the near-field region. With our designed narrow-slit depth, slit dielectric, 
and inter-slit distance, we show that one of the two observation slits can be in a passing 
state while the other one, at the same time, is in a blocking state for a particular wavelength. 
Moreover, at double that wavelength, the previously-blocked observation slit switches to a 
passing state. Therefore, surface-plasmon mediated light-beam splitting becomes wavelength 
sensitive in our model, and can be spatially separated in the transverse direction parallel 
to the array. As a result, it provides us with a unique way for direct color reading in the 
near-field region based on a non-spectroscopic technique. In this paper, we have only cited 
the most relevant and the most recent advances in the fields of light scattering and surface 
plasmons, including works reported by us and other groups. The readers who are interested 
in details of this field are referred to the review article by Garcia- Vidal, et al.^ 

The rest of the paper is organized as follows. In Sec.|2| we derive the scattering- wave 
theory for a non-perfect electric conductor to include the loss of a metallic film in the 
optical-frequency range by employing a surface impedance boundary condition. This theory 
is then applied to study the transmission of an electromagnetic field through an arbitrary 
array of slits perforated on the metallic film and filled with various dielectric materials. 
The issue about using the surface impedance boundary condition for a film was extensively 
discussed in Ref.^. In general, the surface impedance boundary condition can be expressed 
as a linear relation using a (2 x 2) impedance matrix and takes a nonlocal or an integral 
form. The zeroth-order term of the impedance matrix is a local matrix. If the skin depth of 
a metal film is much smaller than its thickness, which is the situation to be considered in this 
paper, the off-diagonal elements of the local impedance matrix, which couple two surfaces 
of a metal film, can be neglected. In Sec. [3} numerical results are presented to demonstrate 
both the passing and blocking states of two observation slits at a particular wavelength and 
at double that wavelength, along with detailed explanations of these two complementary 
states based on the intra-slit and inter-slit dual-wave interferences. The conclusions drawn 
from these results are briefly summarized in Sec.|4j 
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2. Scattering- Wave Theory 

A. Basic Formalism 

In this paper, we consider only p polarization for an electromagnetic (EM) field, written 
as H = (0, Hy, 0), so that the surface-plasmon (SP) wave on a metal-air interface can be 
excited. Here, the H field is assumed translationally invariant in the y direction for the array 
of slits, shown in Fig.[T] We further denote the scalar magnetic-field amplitude as u{x, z) = 
Hy{x, z). The corresponding electric field can be calculated from E = ^/(cjeoes) V x H, 
where uj is the angular frequency and es (real) represents the relative dielectric constants 
of the host materials on the left- (s = L) and right-hand (s = R) side of a metal film. 
Since the metal film, which contains a finite slit array, will be treated as a non-perfect 
electric conductor (non-PEC), we need to employ the so-called surface impedance boundary 
conditionlMSQ] (gmC) for the total EM field. In our current model, the SIBC requires 
du{x, z)/dx = rjsu{x, z) on a metal surface, where rjs = ± /coCs/I^a/— eM(i^) ], ko = uj/c is 
the wave number of the EM field in vacuum, eM(i^) (complex with optical loss) is the metal- 
film dielectric function, and the ± signs refer to the left (minus) and right (positive) surfaces 
of the metal film shown in Fig.jTj The frequency-dependent dielectric function eM(i^) for the 
gold film is obtained by interpolation from the data in the paper by Johnson and Christy'^. 



By applying the SIBC on the left-hand side (L) of the slit array {x < —d), we get 



du{x, z) 
dx 



riLu{x, z)\x=-d-o left surface zj + ij < z < zj^i 



-d-O 



!££, du{x, Z) I 

Ki dx 'x= 



'd+0 



middle slit \z — Zj\ < ij 



(1) 



, ±A^ is the slit index, Zj and 



r]Lu{x, z)|^.=_rf_o right surface Zj^i + < z < Zj — ij 

where 2d is the thickness of the metal film, j = 0, ±1, ±2, 
2ij are the center position and the width of the jth slit, and (real or complex) is the 
dielectric constant of the material inside the jth slit. Similarly, after applying the SIBC to 
the right-hand side (i?) of the slit array {x > d), we acquire 



du{x, z) 



dx 



x=d+0 



rijiu{x, z)\x=d+Q left surface zj + ij < z < zj+i — ij^i 
middle slit \z 



tji du(x, z) 

Kj 



Zj\ < Ij 



dx \x=d—0 

rifiu{x, z)\x=d+o right surface Zj^i + < z < Zj — ij 



(2) 



If we set ?7l = r^R = in Eqs. ([T]) and we will simply get the corresponding boundary 
conditions for a PEC'^. In addition, the continuity of u{x, z) is needed for each slit en- 
try: i.e., 2;)L=_rf_o = u{x, 2)L=_rf+o and z\^^_q = u{x, z)\^^^^q. For Lamellar 
metallic gratings, the PEC boundary condition for the slit side walls was used in calcula- 
tions of plasmon-mediated light reflection, and the results agreed very well with experimental 
result^^'^. Therefore, for the PEC the field normal derivative must be zero along the slit 
side walls, i.e.. 



du{x, z) 



dz 



for all slits \x\ < d . (3) 



In Region I (the left-hand side of the slit array), the total field, including both the incident 
and reflected ones, can be written as'221 



u^'\x, z) = heL I ^ [Gs{P)cos{Pz) + zGa{(3)sm{Pz)] e'^-^P^^^+'^^ 



oo 

- k,eL I ^ [A,{/3) cos{Pz)+zAM sm(/3z)] e-^'=^(^)(^+'^) , (4) 



where Gs{P) and Ga{P) are the symmetric and anti-symmetric spectral components of an 
incident Gaussian beam, given by 



GsiP) = [GniP) + GM] cos(/3^g) - ^ [GnW) - GM] MP^g) , 
GaW) = [GnW) - GpiP)] cos(/3^g) - ^ [Gn{P) + GM] MP^g) , 



(5) 



zg dennotes the Gaussian beam center position, and Gp{(3) and Gni(3) in Eq. (|5| are defined 



as 



exp 



^7'(/3 + /3o 



\2n 



g'iP-Po) 



{niko - 



(6) 
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In Eq. ([6]), Q{(3) is the unit step function, (3o = nikosmOQ, Oq is the incident angle of the 
beam, ul = and ki{j3) = \/'n\}i^ — jS'^ can be either real or complex with lm[A;i(/3)] > 0. 



In a similar way, we find that in Region III (the right-hand side of the slit array) the 
transmitted field takes the forml^^l 



z) = hen j ^ cos(/3^) + sin(/3^)] e^^^(«(-'^) , (7) 





where hr = y/en and k2{/3) = \/n\kl — f}"^ with lm[A;2(/3)] > 0. 



Finally, in Region II (middle slit array), using the eigenmode expansion [subjected to the 
boundary condition in Eq. (Isl)], we obtain '221 



j n '^^^'^ 

X cos[^L(^ + k„ e-^"(^+^) - 6i„ e--^^^-'^)] sin[C(^ - %)] | , (8) 

(Tin ^ -I J 

where n = 1, 2, ■ ■ ■ is the eigenmode index, = (ir/ij) {n — 1) and ^{^ = (vr/£j) (n — 
1/2) are for symmetric and anti-symmetric slit eigenmodes, respectively, and (T^„ ^„ = 



^j^o " {^sn,an)'^ cau be either real or complex, with Re[cr^„ ^j^] > 0. 

When n = 0, the lowest symmetric eigenmode in Eq. ^ corresponds to a uniform EM 
field distribution in the z direction within each slit. The evanescent waves can exist only 
when ki{f3) or k2{f3) is purely imaginary, i.e., (3 > riLko for the reflection side or (3 > riRko 
for the transmission side. On the other hand, the condition for a pure scattered surface 
wave is obtained through ki{P) — ?■ or k2{P) — )■ 0. This condition can be met by /3 = riLko 
for reflection or /3 = nnko for transmission. 
B. SIBC Constraints 

By using the derivative boundary conditions in Eqs. ([T]) and ([2]) at x = ± c/ and the orthogo- 
nality of the continuous Fourier expansions in Eqs. Q and ([T]), a set of constraint equations 
for the unknown Fourier coefficients As{f3), Aa{f3), Bs{f3) and Ba{f3) can be obtained. 

At X = —d, from Eq. ([T]) we get 
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J d(3 {[Gs{(3) + M(3)] cos(/3z) + i [GM + AM] sm{(3z)} 



oo 

-tVLeL J ^ {[GsW) - Am cos(/3z) + t - A^] sin(/3z)} (9) 



for z values within the non-sht regions of a non-PEC, and 



{[GsiP) + A,m cos(/3z) + t [GaW) + AM] sm{Pz)} 
= J2 m - k - ^.1) E + ^sne'"-'') cos[eL(^ - Zj)] 

j n 

+t (aL + e''-'-'') sin[C(^ - z,)]] (10) 
for the sht regions. Similarly, at a; = ci we find from Eq. ^ 



dl3 [BsilS) cos{l3z)+iBail3) sin(/3z)] 

oo 

= -^VReRj ^ [Bs{/3) cos{/3z) + iBa{/3) sm{/3z)] (11) 



for z values within the non-slit regions, and 

oo 

j d/3 [BM cos(/3z) + zBM sm(/3z)] = ^ Q{ij - \z - z,\) [{Hn 

^ 

+aL e^'^-'') cos[6„(^ - z,)] + i [^n + <n e"-'-'") sin[C(^ - z,)]] (12) 

for the slit regions. 

C. Projection of SIBC 

Since the combination of Eqs. ^ and (10) extends over the left surface of a metal film, we 



can project out the symmetric and anti-symmetric Fourier coefficients in u^^\x^ z) through 
multiplying these two equations by cos(/3'z) or sm.{P' z) and integrating over z afterwards. 
This leads to 
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^(/3) + a(/3) = E(E7 [(' 

n { j 



<5an(/3) sin(/3zj 





oo 



- [Pe(/3, /S') + W,i/3, /3')] [Ga{/3') - AaiP')] , 



(13) 



+ f + e^^'^^^"'^) QL(/3) cos(/3^, 



t a 



"an ' "an 
oo 



m / [^a(/3, /?') + H^a(/3, /?')] [Gal/?') " Aa(/3')] 



where the definitions of QinW) and QiniP) can be found from Appendix A. 



(14) 



For the same reason, using Eqs. (11) and (12) we can also project out the symmetric and 



anti-symmetric Fourier coefficients in u through multiplying them by cos(/3'z) or 

sin{P'z) and doing a follow-up z integration. This yields 



Mm = E I [( 

n L j 



bi^ + ai^e'^^i^A Qi^{f3) cos(/3z,) 



QUP) MPz,)] } 



d(3' 





oo 



HP') 



+ Vr 



HP') 



- [P,(/3, /3') + iy,(/3, /?')] i?a(/3') 



(15) 
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+ 




2iai„d 



- [Pe(/3', /3) + W^c(/3', /?)] 



- ^VR 



d/3' 



k2W 



(16) 



The definition of the couphng matrices P,(/3, /?'), ^a(/3, /?')> ^c(/3, P'), W.s{P, /?'), W^a(/3, 
and Wc{P, P') can be found in Appendix B. Without loss of generahty, we assume an order 
for the sht array — oo < (-Z-at — ^-n) < (-S-at + ^-a?) < ■ • ■ < {z_i — < {z_i + < 
(2:0 - io) < {zo + io) < {zi - £1) < (zi + ii) < ■■■ < {zn - iN) < {zn + ^n) < 00. 

\{ TjL = TjR = Q for a PEC, we can exphcitly express!^ the continuous Fourier expansion 
coefficients ^a(/3), Bs{f3) and Ba{f3) by the discrete Fourier expansion coefficients 

a 



bl^, and through Eqs. (13), (14), (15) and (16). If there exists only one slit 



{j = and zq = 0), the PEC single slit does not couple symmetric modes to anti-symmetric 
ones^. In addition, we find Ws{(3, (3') = Wa{/3, /3') = Wc{(3, (3') = and Pc{/3, (3') = 
in this case. Therefore, the non-PEC single slit cannot couple symmetric modes to anti- 
symmetric ones. For a symmetric distribution of slits with respect to Zq = 0, we always 
have — ^ — Ar| = + i^, leading to Pc{f3, /3') = 0. However, in this multi-slit case, 
Wc{P, P') 7^ 0. As a result, symmetric and anti-symmetric modes are coupled to each other. 



If we treat the SIBC in Eqs. (13), (14), (15) and (16) as a perturbation for small r]L and r]fi, 



we can express As{P), Aa{P), Bs{P) and Ba{P) by the discrete Fourier expansion coefficients 
a^„, and to the leading order of the perturbation. Alternatively, if only the 



diagonal contributions for the term proportional to rfi or rf^ are 

keptEnfflS] in Eqs. (pi) 



(14), (15) and (16) [i.e., including only the terms with /?' = ±/3 in the integrals], y4<j(/3), 
Bs{P) and Ba{P) can also be expressed by the discrete Fourier expansion coefficients 
a^„, a{„ and However, such a simplification'^inilSl needs to be justified physically. 



The projected SBIC in Eqs. (13)-(16) leads to the approximate field equations derived in 



Reference j^ * ^^ * ^^ * after applying the "diagonal approximation" for simplification. It is clear 



from Eqs. (13)-(16) that such a "diagonal approximation" can only be justified when the 



coefficients Ps{f3.f3'), Pc{P.f3'), Pa{f3.f3'), as well as the coefficients Wd^.^'), 
Wa{P-P'), are either peaked around /3 = /3' or negligibly small. We further realize from 
Appendices B & C that Ps(/3. /?'), Pa{(3- /?') always peak around (3 = (3' due to the existence 
of the — /?') terms; Ws{P. /?'), Wa{(3. f3') can peak around (3 = (3' only if neighboring slits 
are well separated from each other in comparison with the inverses of /3 and /3' due to the 
existence of sinc-function terms; both Pc{f3. (3') and Wc{f3. f3') can be neglected only for very 
large (3, (3' and (3 ^ (3' due to the existence of cosine terms. Therefore, the advantage of the 
current theory is its immunization from such restrictions and it can be applied to general 
cases. 

D. Integral Equations 

In order to get the closed-form integral equations for a non-PEC, we look to express a{^, 
ai^, V,^ and 6^ in Eqs. (fll, (111, ^ and ^ by and This 



can be achieved by using the orthogonality of the discrete Fourier expansion in Eq. ([8j) and 
multiplying both sides of the continuity conditions for slit entries and exits, i.e. u^^\x = 
—d, z) = u^-^^\x = —d, z) and u'^™^{x = d, z) = u^^^\x = d, z), by cos[^i^,{z — zji)] or 
sm[C,l^,{z — Zj/)], which is followed by an integration of z over all the slit regions. This yields 



J \ sn "sn 



d(3 

km 



+ i 



d(3 



[GM - AM] QW) sm(/3z 



(17) 



7 \ "'an ^an ^ 

oo 

df3 



d(3 



km 



(18) 



oo 



+ i 



d(3 

km 



Bai(3)Qim sin(/3z 



JJ 5 



(19) 
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- V 



d(3 



+ 



dl3 



(20) 



where Xn = 2 for n = 1 and Xn = 1 for n ^ 1. Equations (17) and (18) connect the jth-sht 
field at the entry edge to the total field in Region I, where the forward- moving amplitude, 
or al^, and backward-moving amplitude, or b{^, can be viewed as two independent 
interfering waves with phase delays of 2iai^d or 2ia{^d. The same arguments can be applied 



to Eqs. (19) and (20) at the exit edge of the jth slit. 



Equations (17) through Eq. (20) can be formally solved analytically, which leads to 



y (2) _ ^-2iai„d y (1) 



2i sin(2o"snrf) 
2i sm{2aind) 



V 



yW 

jn 



2i sm{2aind) 

^ ^^jn ^^jn 

2i sm{2aind) 



(21) 



Here, and represents the forward-moving waves, while and represents the 
backward-moving waves. Whenever 2ai^d/7T or 2ai^d/7c becomes an integer, the constructive 
or destructive dual-wave interference will occur at two the edges of the jth slit for the nth 
eigenmode. In this case, the slit behaves either like a passing filter or like a cavity for field 



trapping. The definitions of yjnK ^jni -^"jn ■> ^^"^ -^)n introduced in Eq. (21) can be found 



'(1) 



(2) 



in Appendix C. Here, X^^^ and Y.-^^^ come from the contribution of u^'^\x, z) at the left slit 



^(1) 



entry, while X-^^ and Y^'^'' come from the contribution of ^'•^^^^(x, z) at the right slit exit. For 



'(2) 



a PEC, by substituting Eqs. (|13j), (|14j), (|15j) and (jiej) with r]L = r]R = Q into Eqs. (|17j)-(l20j), 
we get a set of inhomogeneous linear equations'221 with respect to a^„, 6^^, and 



Using the results in Eq. (21), from Eqs. (13), (14), (15) and (16) we finally obtain four 
integral equations for As(/3), Aa{P), Bs{f3) and Ba{f3) 
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oo 


oo 

+ riLj [Pc(/3, n + P^c(/3, /?')] - Mn 



oo 





oo 





oo 





oo 





oo 

GaW) + Aa{(3) = -rjL j [PcW, /3) + W,{(3', (3)] [Gs{P') - 



oo 





oo 

/■ d/3' 





oo 





oo 





oo 



km 





12 



oo 



+ ^7^ 





oo 

d(3' 



^ [P,(/3, /?') + W,{(3, (3')] 



oo 



oo 





oo 

dp' 



HP') 



oo 





+ iMi'\p,p') [Gs{P')-As{P')]] , (24) 

oo 

BaW) = -^i^ 1 ^ [Pc(/3', /?) + W^c(^', /?)] 


oo 

- ^VR J ^ [PaiP, P') + W^a(/5, P')] B^iP') 





oo 

dp' 



HP') 



oo 

dp' 



- [Ni'\p, P') B,(P') + iN('\p, P') B^iP')] 





oo 





oo 



+ tMi'\p,p') [Gs{P')-A,{P')]} . (25) 
The expressions for the matrices ri°'(^, P'), Ri°'\p, P'), M^°'\p, P') and#(^, P'), as well 
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as for Tr'{l3, /3'), Rr{/3, /3'), Mr{f3, f3') and Nr{f3, f3'), with a = 1, 2, 3, 4 are given ex- 
plicitly in Appendix D. In general, /?'), Rs':a\^, ^ft'^/?, and N^s!a^\f3, 
come from the contribution of z), while Ti,^a^^(/3, P'), M^%^\p, and 
come from the contribution of u^^\x, z). If the array is infinite and all the slits, 
as well as the filled dielectric materials, are identical for a periodic system, P^'^ * ^" * ^^ * ^^ ', the con- 
tinuous variable /3 becomes a discrete reciprocal wave number j{27T/V) where T> is the array 
period and j = 0, ±1, ■ • • . 

There are three different interferences discussed in this paper, namely, intra-slit and inter- 
slit dual- wave interferences and Fabry-Perot sole- wave interference. First, for the intra-slit 
dual-wave interference, we consider the interference between the forward and backward 
moving slit waves. For the backward moving wave, there already exists a vr phase delay 
and the latter interferes with the former after its reflection from the entrance side of a slit. 
Therefore, the constructive interference condition (on-state) between these two waves is 
simply given by ^/KkoAd = (2m — l)7r, where ko is the wave number in vacuum, n is the slit 
dielectric constant and m = ±1, ±2, ■ ■ ■ is an integer.For the same reason, the destructive 
interference condition (off-state) is given by ^/KkQ4d = 2rmT. Second, by the inter-slit dual- 
wave interference we mean the interference between a backward moving reflected slit wave 
and a surface wave propagated from another slit and entering forward into the side slit 
studied. In comparison with the intra-slit dual-wave interference, there exists an additional 
phase compensation for the forward moving slit wave in this case, which depends only on 
the slit separation for the surface wave but not on the metal dielectric function. Therefore, 
we expect a strong influence from the inter-slit dual-wave interference on the intra-slit dual- 
wave interference. At last, for the Fabry-Perot sole-wave interference, on the other hand, 
it is the interference between a forward moving slit wave and the same wave after it has 
been reflected twice successively by the exit and entrance sides of a slit. Therefore, the 
constructive Fabry-Perot interference condition takes the form of ^/K,kQAd = 2rmT, which 
sits at the same position as that of the destructive intra-slit dual-wave interference, and its 
strength goes up with increased finesse in a longer slit. 

3. Numerical Results and Discussions 

In this paper, we would like to demonstrate a direct optical reading of fundamental 
and second-harmonic near-field photon emissions in the near-field region using a non- 
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spectroscopic technique with specifically designed slit depth, slit dielectric material, and 
inter-slit distance. In our calculations, we take €1 = 1 (air) and = 20, 25 (high-dielectric- 
constant oxides). For a triple-slit structure, we illuminate the middle slit and find conditions 
under which one of the side slits is in a passing state for a particular wavelength (with the 
other slit in a blocking state) while the other side slit is in a passing state for twice that 
wavelength (with the other slit again in a blocking state). Here, the passing state of a slit 
refers to the fact that a surface wave propagated from the middle slit can pass through 
a vertical slit to get to Region III from Region I, while the blocking state corresponds to 



a surface wave not passing through a slit due to interference effects. From Eqs. (17)-(20) 
we find that the forward and backward moving waves inside a slit can be treated as two 
independent waves, which are decided by the fields in Regions I and III, respectively. As 
a result, the intra-slit interference in this paper can be viewed as a dual-wave interference, 
which is quite different from the geometric series result of a typical Fabry-Perot (FP) slit 
model. FP cavity analysis utilizes intra-slit reflection and transmission coefficients for a 
single wave which are not related to the total fleld structures at the slit edges. It is also 
important to know that there is already an extra vr-phase shift in Eq. (|8| due to opposite 
signs for the forward-moving wave (+) and the backward-moving wave (— ). Therefore, for 
a given incident wavelength Aq in vacuum, the constructive intra-slit dual-wave interference, 
or the passing state of a slit, is found to satisfy the condition: 4d = (2m — 1)Ao/2-\/k with 
m = 0, ±1, ±2, ■ ■ ■ , where k is the slit dielectric constant. Moreover, the condition for 
the destructive intra-slit dual-wave interference, or the blocking state of a slit, is given by 
4d = 2mAo/2A/K with m = ±1, ±2, ■ ■ ■ . Whether for plane wave or SPP excitation of 
a single slit, the constructive or destructive dual- wave interference results agree with the 
zero-order transmission coefficient in a single-mode approximation!^. The above interfer- 
ence conditions are not directly derived from the theory described in Sec. 2. However, these 
phenomenological arguments are found to explain the numerical results in this paper pretty 
well. 

By employing the mechanism for the intra-slit dual-wave interference discussed above. 
Fig. [2] simultaneously displays the passing state of the left slit as well as the blocking state 
of the right slit at Aq = 0.588 fim, where the middle slit is used for a local front-side SP 
excitation, as seen from Eq. (|6]), and is filled with a wide-band attenuator to prevent light 
from leaking into Region HI. When the slit depth 2d becomes larger than Aq and the slit 
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width 2ij is only half of Aq, all the slit modes for p polarization become strongly attenuated 
except for the lowest symmetric one which has a uniform distribution in the transverse {z) 
direction of a slit. For this lowest symmetric mode, the dual- wave constructive- interference 
condition for the left-slit is satisfied due to 4d = 7Ao/2^k_i and = 1- We further find 

that the dual-wave destructive-interference condition for the right-slit is met at the same 
time due to Ad = 14Ao/2y^ and = 2. This fully explains the observed left-slit passing 
state as well as the right-slit blocking state in Fig. [2] 

In order to get a complete picture about the intra-slit dual-wave interference after the 
propagation of an SP wave locally excited at the middle slit, we present in Fig. |3] the aver- 
aged ratio of \Hy{x, z)\'^ over a slit at Aq = 0.588 /im as a function of slit depth d for two 
observation slits at a distance 6 slightly away from the backside {z = d) of the metal film. 
For the left slit at z = Z-i (blue solid curve), we find from Fig. |3] that there exist four minima 
at d/Xo = 2/8, 4/8, 6/8 and 8/8, which agree with the dual-wave destructive-interference 
condition, i.e., 4d = 2mAo/2^/t_i with m = 1, 2, 3, 4 and = 1. Similarly, for the right 

slit ai z = Zi (red dashed curve), we find eight minima at d/Xo = 2/16, 4/16, ■ ■ ■ , 14/16 and 
16/16, which also agree with the right-slit dual- wave destructive-interference condition, i.e., 
Ad = 2mXo/2y/Ki with m = 1, 2, ■ ■ ■ , 7, 8 and ^/ki = 2. In addition, a maximum is always 
seen between two adjacent minima for the right slit, which meets the dual- wave constructive- 
interference relation, i.e.. Ad = (2m — l)Ao/2y^ with m = 1, 2, ■ ■ ■ , 7, 8. However, for the 
left slit, instead of four maxima as expected, we only see four dips at d/Xo = 1/8, 3/8, 5/8 
and 7/8 due to the effect of an inter-slit dual- wave interference explained below. We note 
that the slit widths in Figs.|2]and|3]are different, but it will not change the peak and valley 
positions for the dominant lowest symmetric slit mode. 

It is important to mention that the inter-slit distance in Fig.[3]is set to be (zi — z_i)/Ao = 
3.75, which is a quarter wavelength smaller than a multiple of the wavelength Aq. In order 
to understand the inter-slit dual-wave interference effect, we show in Fig. |4] the slit-averaged 
ratio of \Hy{x, z)\'^ at Aq = 0.588 /im as a function of d for z_i = — 1.5Ao in Fig.|4|^a) and 
z^i = -2Ao in Fig.gb) with zi = 2Ao fixed. When d/Xo = 1/8, 3/8, 5/8 and 7/8, from 
Fig.[3]we already know that the contributions from the exit-edge reflected forward-moving 
waves are out of phase at the left and right slits since the former is in a passing state while 
the latter is in a blocking state. Therefore, for the case with {zi — z^i)/Xo = 3.5 in Fig.|4]^a), 
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we expect a constructive inter-slit dual-wave interference to occur at the left slit (blue solid 
curve), which can be verified by noticing the change of a dip in Fig.|3]into a maximum in 
Fig.|4]^a). However, when {zi — z^i)/\q = 4 as shown in Fig.|4](b), the dip in Fig. |3] changes 
into a full minimum (blue solid curve) as a consequence of the destructive inter-slit dual-wave 
interference. 

As a complementary result to the left-slit passing state at a particular wavelength Aq = 
0.58 fim in Fig.|2| we present another contour plot for the EM field distribution at double 
that wavelength Aq = 1.176 /im in Fig.|5j where the right-slit is in the passing state in this 
case. Here, the constructive intra-slit dual-wave interference condition Ad = 7Ao/2^//€7 has 
been met for the right slit. However, the left slit is neither in the passing state nor in the 
blocking state at Aq = 1.176 fim, which is seen in an intermediate state between these two 
with Ad = (7/2)Ao/2^/KrT. 

To get a complete picture about the intra-slit dual-wave interference at Aq = 1.176 /im, 
we show the slit-averaged ratio of \Hy{x, z)\'^ at Aq = 1.176 yum in Fig. [6] as a function of 
d. Indeed, we find from this figure that there exists a set of minima for the right slit at 
2d/Xo = 2/8, 4/8, 6/8 and 8/8 with ^/k^ = 2, as well as a group of maxima at 2d/Xo = 1/8, 
3/8, 5/8 and 7/8. We also notice here that dips, instead of maxima, still show up for the 
left slit at 2c// Ao = 2/8 and 6/8 with = 1 due to incomplete inter-slit dual- wave 

constructive interference as discussed in connection with Fig.|4j Interestingly, we also find 
very narrow peaks right above two of the four blocking states of the right slit, and these 
sharp peaks get stronger with increased slit depth. The occurrence of two extremely sharp 
peaks in this figure can be qualitatively attributed to the result of the Fabry-Perot sole- 
wave interference between the entry and exit edges of the right slit with a finite value for 
the average reflection coefficient (or a large finesse). The peak positions from a constructive 
Fabry-Perot interference coincidentally overlaps with two of the four slit blocking states. 
Although the slit widths in Figs.[5]and|6]are different, it will not change the peak and valley 
positions for the dominant lowest symmetric slit mode. 

The current paper deals with a real metallic film by imposing the SIBC in Eqs. Q and 
(|2| for a finite metal conductivity. This facilitates the propagation of the SP polariton wave, 
which is excited at the middle slit, to two neighboring side slits. Consequently, for a PEC 
with rjL = rjn = we expect the passing state of the left slit in Fig. [2] will become much 
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weaker, which can be clearly seen from the calculated transmitted near-field distribution in 
Fig.[7| In the case of a PEC, no SP polariton wave is excited on the backside of the metal 
film, and the refiection from the illuminated middle slit on the front-side of the film becomes 
very collimated in the near-field region. Moreover, the maximum intensity in Fig.[2] occurs 
farther from the film than it does in Fig.[7| and the angular distribution of the intensity of 
the transmitted field is broader in Fig. |2] than it is in Fig.[7j There exists a major difference 
in the coupling between slits on the interface of a PEC or a real metal. For the former, 
only radiative modes out of the metal plane can contribute, while both surface-plasmon and 
radiative modes will contribute to the latter. 

Finally, we know from Fig.|5] that the left slit is in an intermediate state although the 
right slit is in a passing state at Aq = 1.176 fim. The ideal situation is that the left slit could 
be forced into a blocking state at this wavelength. This goal can be reached if the filled 
dielectric medium in the left slit can be tuned from y^/tTY = 1 to ^/k-i = 4, which has been 
simulated by the calculated transmitted near-field distribution displayed in Fig.[8j 

4. Conclusions 

In conclusion, motivated by the previous Green's function formalism'^^ini^ we have derived 
a scattering-wave theory to study the wavelength-dependent detection of surface-plasmon 
mediated light-beam splitting into two by a triple-slit structure perforated by a gold film and 
filled with different slit materials. This can be viewed as a new addition to the previously 
demonstrated longitudinal color-dependent light focusing using a finite groove array with 
various groove widths in a parabolic pattern. For a specifically chosen slit depth, filled slits 
dielectric material, and inter-slit distance in the deep sub-wavelength regime, we have found 
that only one of the two side observation slits is in a passing state for a particular wavelength, 
but the other blocked slit switches to a passing state at double that wavelength. In this sense, 
surface-plasmon mediated light-beam splitting becomes wavelength sensitive, and a single 
light-beam incidence with two wavelengths can be separated along the transverse direction 
parallel to the array. This provides us with a direct optical reading in the near-field region 
based on a non-spectroscopic technique. 
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Appendix A 



The overlap integrals initially introduced in Eqs. (13) and (14) are defined as follows 



dz cos{C,l^z) cos(/3z) 



sinc[(/3-ej^.]+sincp + ej^. 



(26) 



dz sm{^i^z) sin(/3z) 



smc 



a)£,]-sinc[(/3 + eL)^,] 



(27) 



and sinc(x) = sinx/x. 



Appendix B 



The coupling matrices initially introduced in Eqs. (13), (14), (15) and (16) are defined as 



— k-jv— ^-jvl oo 

— oo ZN+iN 

7r[5(/3-/3') + 5(/3 + /3')] 

\Z-N - ^-Af| 



Zn + iN 



{sinc[(/3 + /3')k-7v - i-N\] + sinc[(/3 - /3')\z.n - ^-7v|]} 
{sinc[(/3 + /3'){zN + In)] + sinc[(/3 - /3')(^^ + ^n)]} , 



(28) 
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— \z-n—1-n\ oo 

Pa(^, p') = j sm{^'z) sm{pz) dz + J sm{j3'z) sm{pz) dz 

-oo ZN+tN 

+ {smc[(/3 + /3')k-iv-^-iv|]-smc[(/3-/3')k-iv-^-iv|]} 

+ ^^^^ {smc[(/3 + P'){z^ + e^)] - smc[(^ - P'){zr, + i^)]} , 



— \z-n—£-n\ oo 

P^(/3, = y" sm(^'^) cos(^z) + J sm{p'z) cos(^z) dz 

— oo zm+<^n 

1 rcosp + /3')k-iv-^-iv|] cos[(/3-/3')k-iv-^-iv| 



+ 



2 { (3 + (3' (3 -(3' 

1 / cos[(/3 + /3')(-2iv + ^iv)] cos[(/3 - /3')(;^iv + ^at)] 



(/?, /?') = XI / cos(^'z) cos(/3z) dz 



2 



+ sinc[(/3' + /3)(z,+i-£,+i)]) 

- (z,. + £j) (sinc[(/3' + /3)(z, + £,)] + sinc[(/3' - /3)(z,. + £,)])} , 



7V-1 ^i+i-^i+i 

Wa{(3, f3')= J2 J ^'^i^'^) ^^i^^) 

1 

+ (z, + i,) (smc[(/3' + /3)(z, + - smc[(/5' - /3)(z, + i,)])} , 
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N-1 -^+^ 

WciP, P')= Yl / sm(/3'2) cos{/3z) dz 

^ AT 'J 



j=~N 



j=-N 

cos[(/3' - fi){z, + £,)] ^ cos[(/3' + + ^,)] 



cos[(/3' - - cos[(/3' + - 



Appendix C 



We have introduced in Eq. (21) the following amplitudes 





X {[Gs{l3) - A,{I3))] cos(/3^,) + i - sin(/3^,-)} , 



= — [l^, QlnW) [BsW) cos(/3z,) + iBM sm(/3^,)] 









X - A^m cos(/3z,) + 2 [G,{I3) - sin(/3^,)} 



K2) (^in^R f d(3 







Appendix D 



In Eqs. (22) through (25) we have introduced the following matrices 



21 



X 



X 



X 



R^a\/3, /?') 
R^a\P, P') 
R^a\f3, n 



-E 



XnK'j sm{2aind) 



en cos{l3zj) cos{l3'zj) 
cos{P Zj) sm{P' Zj) 
eLC0s{2ai^d) cos{/3zj) cos{f3'zj) 
ei cos{2ai^d) cos{(3zj) sm{(3' zj) 



ITT 



£■ 



Xni^j sm{2aind) 



en sni{p Zj) cos{P' zj) 
€r sm{(3zj) sm{f3'zj) 
ex, cos(2(T^„(i) sm{l3zj) cos{(3'zj) 
cos(2(j^^(i) sm{l3zj) sm{(3'zj) 



er sva{l3zj) cos{f3'zj) 
tR sm{f3zj) sm{f3'zj) 
El cos(2cr^„(i) sm{f3zj) cos{f3'zj] 
eLC0s{2ai^d) sm{(3zj) sm{(3' Zj) 



Xn^^i sin(2cri„(i) 
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X 



€r C0s{/3Zj) C0S{P' Zj) 

er cos{ f3zj) sm{l3'zj) 
ei, cos(2o'^„(i) cos{f3zj) cos{f3'zj] 
eLC0s{2al^d) cos{(3 Zj) sm{P' Zj) j 



ivP(/3, /?') 
ivP(/3, /?') 
N^iP, (5') 



n,3 



Xnf^j sm{2aind) 



QLW)QLW') 



X 



ejiCos(2ai^d) cos{(3zj) cos{f3'zj) 
cos(2(T^„(i) cos{f3zj) sm{f3'zj) 
El cos{(3zj) cos{(3' Zj) 
€l cos{/3zj) sm{l3'zj) 



mP{/3, (3') 



n,3 



Xni^j sin(2ai„ci) 



QM QLW) 



X 



eRCOs{2al^d) sm{f3zj) cos{f3'zj) 
€RCos{2al^d) sm{/3zj)sm{l3'zj) 
€l sm{/3zj) cos{/3'zj) 
€l sm{Pzj) sm{f3'zj) 



-E 

ITT 



sn^J 



XnK'j sm{2aind) 



QinW)QlnW') 



X 



COS (2(7^^0?) sm{f5zj) cos{(3' Zj) 
6RCOs{2ai^d) sm{Pzj) sm{P' Zj) 
€l sm{f3zj) cos{P' Zj) 
El sin(/3^j) sva.{l3'zj) 
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MPif3, (3') 



ITT 



£■ 



XnK'j sm{2aind) 



QiniP) QiniP') 



X 



€RC0s{2ai^d) cos{/3 Zj) cos{l3' Zj) 
eRC0s(2a{^d) cos{p z j) s\n{p' Zj) 
€l cos{/3zj) cos{/3'zj) 
El cos{(3zj) sva.{l3' Zj) 



(45) 
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Fig. 1. (Color online) Illustration for a z direction slit array (brown) which extends in the 
y direction, where Zj and 2ij are the center position and the width of the jth slit with 
j = 0, ±1, ■ ■ ■ , ±A^. The regions at the left- and right-hand side of the slits are denoted as 
Region I and Region III, respectively, with real dielectric constants ei and e^. The region 
for the slit array is denoted as Region II, and slits are filled with medium having a dielectric 
constant Kj (real or complex) for j = 0, ±1, ■ ■ ■ , ±A^. The depth of slits in the x direction 
is 2d, and €4/(0;) represents the dielectric function of the metal film containing slits. A 
Gaussian beam is incident on the slit array from the left side with an incident angle 6q and 
at a center position z = zq. The incident wave number is yfeik^ and /3o is the incident wave 
vector along the z direction. 
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z(|Liin) 

Fig. 2. (Color online) Contour plot of \Hy{x, z)\'^ for p-polarization normal incidence (from 
upper surface) with 9q = 0°, where ra^ = 1 and nji = 4.5. In our calculations, we set the 
parameters as follows: £_i = io = ii = C/4, = 1, ^/Eq = 1 + 30i, = 2, z_i/C = 

— 1.75, zq = zg = 0, zi/( = 2, d = (7/8) C, 9 = G^o, and Aq = where ( = 0.588 /im. 
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Fig. 3. (Color online) Plot for calculated Tj = (l/2£j) / dz \Hy{d + 5. 

(blue solid curve) and 1 (red dashed curve) as a function of slit depth d/Q for p-polarization 
normal incidence as in Fig.[2| where ul = 1, ur = 4.5, S = C/2 and the vertical black dashed 
lines indicate the positions determined by d/( = 2/8, 4/8, 6/8, 8/8. In our calculations, we 
set the parameters as follows: £_i = io = h = C/6, ^/t-i = 1, a/kq = 1 + 30i, = 2, 
z_i/( = —1.75, zq = zg = 0, zi/( = 2, g = 6io, and Aq = (, where ( = 0.588 /im. 
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Fig. 4. (Color online) Plot for calculated Tj = {l/2ij) J dz \Hy{d + 6, z - Zj)\'^ for j = -1 

(blue solid curve) and 1 (red dashed curve) as a function of slit depth d/( for p-polarization 
normal incidence as in Fig.[2| where = 1, = 4.5, 5 = C/2. In our calculations, we 
set the parameters as follows: £„i = io = ii = C/6, = 1, = 1 + 30i, = 2, 

-2^0 = -2^G = 0, zi/( = 2, g = 6io, and Aq = where ( = 0.588 /im. Here, we chose 
z^i/( = —1.5 [in (a)] and z_i/C = —2 [in (b)], respectively. 
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Fig. 5. (Color online) Contour plot of \Hy{x^ z)\'^ for p-polarization normal incidence (from 
upper surface) with 6q = 0°, where ul = 1 and nji = 4.5. In our calculations, we set 
the parameters as follows: £_i = io = ii = C/4, ^/t-i = 1, a/kq = 1 + 30i, = 2, 
z_i/( = —1.75, Zq = zg = 0, Zi/( = 2, d = (7/8) C, 9 = 6£o, and Aq = 2^, where 
C = 0.588 /im. 
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Fig. 6. (Color online) Plot for calculated Tj = {l/2ij) J dz \Hy{d + S, 



Z — Zi 



for j 



-1 



(blue solid curve) and 1 (red dashed curve) as a function of slit depth rf/C for p-polarization 
normal incidence as in Fig.jsj where = 1, nji = 4.5, 5 = C/2 and the vertical black dashed 
lines indicate the positions determined hj d/( = 2/8, 4/8, 6/8, 8/8.. In our calculations, we 
set the parameters as follows: £_i = io = ii = C/6, y/i^-i = 1, ~ 1 + SOi, = 2, 
z-xjC, = —1.75, zq = zq = 0, zi/( = 2, g = 6£o, and Aq = 2^, where ( = 0.588 /^m. 
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z (|im) 

Fig. 7. (Color online) Contour plot of \Hy{x, z)\'^ for p-polarization normal incidence (from 
upper surface) with 6q = 0°, where ul = 1 and = 4.5. In our calculations, we set the 
parameters as follows: = io = ii = C/4, = 1, = 1 + 30i, = 2, Z-i/( = 

— 1.75, zq = zg = 0, = 2, d = (7/8) g = 6£o, and Aq = C; where C = 0.588 /xm. In 
this case, we set ?7l = ?7_r = for a perfect electric conductor. 
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Fig. 8. (Color online) Contour plot of \Hy{x^ z)\'^ for p-polarization normal incidence (from 
upper surface) with 6q = 0°, where ul = 1 and nji = 4.5. In our calculations, we set 
the parameters as follows: = io = ii = C/4, y/^Zi = 4, = 1 + 30i, = 2, 
z-xjC, = —1.75, zq = zg = 0, zi/( = 2, d = (7/8) C, 9 = 6£o, and Aq = 2^, where 
C = 0.588 /im. 
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